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Abstract. We study the Bogoliubov-Dirac-Fock model introduced by Chaix 
and Iracane (J. Phys. B., 22, 3791—3814, 1989) which is a mean-field theory 
deduced from no-photon QED. The associated functional is bounded from be- 
low. In the presence of an external field, a minimizer, if it exists, is interpreted 
as the polarized vacuum and it solves a self-consistent equation. 

In a recent paper, we proved the convergence of the iterative fixed-point 
scheme naturally associated with this equation to a global minimizer of the 
BDF functional, under some restrictive conditions on the external potential, 
the ultraviolet cut-off A and the bare fine structure constant a. In the present 
work, we improve this result by showing the existence of the minimizer by 
a variational method, for any cut-off A and without any constraint on the 
external field. 

We also study the behaviour of the minimizer as A goes to infinity and show 
that the theory is "nullified" in that limit, as predicted first by Landau: the 
vacuum totally cancels the external potential. Therefore the limit case of an 
infinite cut-off makes no sense both from a physical and mathematical point of 
view. 

Finally, we perform a charge and density renormalization scheme applying 
simultaneously to all orders of the fine structure constant a, on a simplified 
model where the exchange term is neglected. 



1. Introduction 

Despite the incredible predictive power of Quantum Electrodynamics (QED) its 
description in terms of perturbation theory restricts its range of applicability. In 
fact a mathematical consistent formulation is still unknown. We want to make a 
tiny step in that direction. 

Following ideas of Chaix and Iracane jx|, we study in this paper a model for the 
polarized vacuum in a Hartree-Fock type approximation. This so-called Bogoliubov- 
Dirac-Fock (BDF) model has been derived from no-photon QED in as a possible 
cure to the fundamental problems associated with standard relativistic quantum 
chemistry calculations. 

The vacuum polarization (VP) is, quoting |16| . "one of the most interesting of 
the phenomena predicted by contemporary quantum electrodynamics" . Although 
it plays a minor role in the calculation of the Lamb-shift for the ordinary hydrogen 
atom (comparing to other electrodynamic phenomena), it is important for High- 
Z atoms [23 0^1 an d even plays a crucial role for muonic atoms I2T] . It also 
explains the production of electron-positron pairs, observed experimentally in heavy 
ions collision P 1321 l3"Tl ITfl . 

In [Jj, Chaix and Iracane noticed that the vacuum polarization effects are "nec- 
essary for the internal consistency of the relativistic mean-field theory and should 
therefore be taken into account in proper self-consistent calculations, independently 
of the magnitude of the physical effects" page 3813]. Taking into account these 



The authors are thankful to Robert Seiringer and Vladimir M. Shabaev for valuable comments. 
They acknowledge support through the European Union's IHP network Analysis & Quantum 
HPRN-CT-2002-00277. E.S. acknowledges support from the Institut Universitaire de France. 



2 



C. HAINZL, M. LEWIN, AND E. SERE 



effects, they restricted the no-photon QED Hamiltonian (normal-ordered with re- 
spect to the free electrons and positrons) to Bogoliubov transformations of the free 
vacuum. This allowed them to obtain a bounded-below energy, a property which 
is a huge advantage compared to the usual Dirac-Fock theory the Dirac-Fock 
energy is unbounded from below, which is the cause of important computational 
El and theoretical Q31 EHH3 EH] problems . 

In this paper, we show the existence of a global minimizer for the Bogoliubov- 
Dirac-Fock functional of Chaix-Iracane in the presence of an external field, which 
is interpreted as the polarized vacuum. This vacuum is represented by a projector 
of infinite rank which solves a self-consistent equation: it is the projector on the 
negative eigenspace of an effective mean-field Dirac operator taking into account the 
vacuum polarization potentials. This equation naturally leads to an iterative fixed- 
point procedure for solving it. In a previous work we proved the convergence 
of such an iterative scheme to a global minimizer of the BDF functional, but under 
some assumptions on the external field and the ultraviolet cut-off. Our goal here is 
to show the existence of a minimizer without any restriction, by means of a direct 
- non constructive - minimization argument. 

In the case where no external field is present, the free vacuum is already known 
to be a minimizer of the BDF energy, as shown by Chaix-Iracane-Lions |S] and 
Bach-Barbaroux-Helffer-Siedentop 0|. In j^, an external field is also considered 
but vacuum polarization is neglected: the model studied there is thus very different 
from the one considered by Chaix-Iracane in |Jj and in the present paper. 

Of course the vacuum case is only a first step in the study of the Chaix-Iracane 
model. In order to consider atoms and molecules, one has to minimize the BDF 
energy in a fixed charge sector, a much more complicated problem from a math- 
ematical point of view. A minimizer would then solve a self-consistent equation 
which takes the form of the usual unprojected Dirac-Fock equations, perturbed by 
the vacuum polarization potentials. 

To deal with divergencies, we impose a ultraviolet momentum cutoff A. Our only 
restriction on A is its finiteness. Additionally we study the behaviour of our solution 
when A — > oo and show that the model becomes meaningless since the vacuum 
density totally cancels the external potential. In physics, this "nullification" of the 
theory as the cut-off diverges has been first predicted by Landau et at |33M32l lTll34 j 
and later thoroughly studied by Pomeranchuk et al. |4(J |. 

We also discuss a simplified model in more detail, neglecting the exchange energy. 
For the corresponding self-consistent solution we perform a fully - to any order in the 
coupling constant a - consistent charge renormalization scheme. This procedure has 
already been performed in perturbation theory by means of Feynman diagrams, see 
e.g. |411 page 194] and In particular, we recover the well-known fact |32U33l lT] 
(see also, e.g., j^Hl Eq. (7.18)]) that the physical (renormalized) coupling constant 
a r is related to the bare a by a relation of the form 



(1) 



1 + aB A 

where B\ ~a^oo 2/(37r) log A. Therefore the limit case of an infinite cut-off appears 
as unphysical |33| since it would correspond to a r = 0, which means no more 
electrostatic interactions. 

The paper is organized as follows. In the next section, we recall the BDF model. 
Our main existence result is stated in Section 3, together with the behaviour of the 
solution as A — > oo. In Section 4, we study the reduced model and interpret the 
self-consistent equation thanks to a renormalization of the charge and the density. 
Finally, the last section is devoted to the proof of our main results. 
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2. The Bogoliubov-Dirac-Fock model 

For the sake of clarity, we first briefly recall the Bogoliubov-Dirac-Fock (BDF) 
model introduced by Chaix-Iracane in [7j and studied in [21] . Details can be found 
in |2S1. 

We use relativistic units h = c = 1, set the particle mass equal to one and 
a = e 2 /(4ir). We emphasize that in the first part e represents the bare charge of 
the electron. We assume the presence of an external field if = n * A describing one 
or more extended nuclei with overall charge density n(x). We do not assume in this 
work that n is a non-negative function, since our model allows to treat the vacuum 
interacting with both matter and antimatter. We denote by D° = a ■ p + (3 the free 
Dirac operator and by D v := D° — atp the Dirac operator with external potential. 
Throughout the paper we use the notation X(-oo,o)(H) to denote the projector on 
the negative spectral subspace of H. In the physical literature X(~oc.o)(H) is often 
denoted as A~(H). 

When the external field is not too strong, a good approximation is to use the 
Furry picture ^H] in the Lamb-shift calculations of atomic bound states (see, e.g., 
[371 135]). This means that, in order to evaluate corrections due to Vacuum Po- 
larization, the dressed vacuum is represented by the projector associated with the 
negative spectrum of the Dirac operator with external potential D v 

In reality, the polarized vacuum modifies the electrostatic field, and the virtual 
electrons react to the corrected field. This remark naturally leads to a self-consistent 
equation for the dressed vacuum of the form 

P sc f = X(-oo-fi)(D' p + V e s) 

where V c s is an effective potential already including the Vacuum Polarization po- 
tentials. The BDF model jj] allows to derive such an effective potential V c s in a 
self- consistent way, P sc f being interpreted as a minimizer in the class of Bogoliubov 
transformations of the free vacuum P° — X(_oo;0)(-O )- 

In practice, V e e can be computed by a fixed point iterative procedure studied in 
details in [2S]. If one starts the procedure from P°, the first iteration gives P v , and 
this explains why the Furry picture is a good approximation. But corrections to 
the Furry picture are necessary for high accuracy computations of electronic levels 
near heavy nuclei. These corrections can be interpreted as the second iteration in 
a Banach fixed-point algorithm (see, e.g., [23 section 8.2]). 

Self-consistent equations leading to a fixed-point iterative scheme are well-known 
and widely used in full QED. The solutions of the Schwinger-Dyson equations |4fil 
112) involving the different four-dimensional Feynman propagators are usually found 
by means of perturbation theory. Our approach for the special case of the Hartree- 
Fock theory without photon studied in this paper is mathematically rigorous, non- 
pcrturbative and works for any charge Z of the external potential. 

The momentum cutoff A is implemented in the Hilbert space 

(2) £ A = {/eL 2 (R 3 ,C 4 ) I su PP /cB(0,A)}, 

that is the space of spin valued functions whose Fourier transform has support inside 
a ball with radius A. Such a sharp cut-off does not allow to keep gauge invariance 
when photons are present. Since we neglect photons, we shall however use (J2J for 
simplicity. 

The space fy\ can be decomposed as a direct sum of the negative and positive 
subspaces of the free Dirac operator D°, i.e. S)a = ® ?)+ where ff_ — P°S ! )a 
and 5° = (1 - P°)F) A , P° = X(-oc;0)(£>°)- The Fock space T is built upon this 
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splitting as usual |71I55]: 

oo 
n,m— 1 

where J 7 ^ := A"=i i s tne "--electron state subspace, ^ m - ) : = AJLi Cfy°_ is the 
m-positron state subspace, and J-^ = J-^ = C. Here C is the charge conjugation 
operator |55j . The bare annihilation operators for electrons ao(/) and positrons 
&o(/) are then defined in the usual way 0EH], for any / S Sja- The field operator 
reads 

*(/)= a (f)+b* (f). 
The free vacuum J7o = 1 6 C C T is caracterized up to a phase by the properties 
ao(/)fio = b (f)Q = for any / <E $)a, and ||Oo||^ = 1- 

Let us now define the BDF class in the Fock space. Given a new (dressed) 
projector P, we define the dressed annihilation operators by ap(f) — ^((1 — P)f) 
and bp(f) = ^f*(Pf). The associated dressed fip is a state in the Fock space 
such that ap(f)ftp = bp(f)£lp — for any / e 5ja, and ||f2p||.F = 1. By the 
Shale-Stinespring Theorem |50|. such an Qp is known to exist and is unique up to 
a phase, if and only if P — P° € 62 (^a), the space of Hilbert-Schmidt operators on 
$)A (see also (HO])- The state fip can be expressed as a rotation of the free vacuum, 
flp = U51o, U being called a Bogoliubov transformation. An explicit formula for 
Q P can be found in lots of papers [SSI EOl HS1 H31 EZl EDI • The BDF class is therefore 
the subset of T 

B := {n P I P orth. projector, P - P° e 6 2 (£a)} ■ 
The charge of Sip can be easily computed 

(3) (n P |G|fi P ) - tr(P (P-P )P )+tr((l-P )(P-P°)(l-P )) 

= tr(Q~)+tr(Q++) 

where Q = P - P° e 6 2 (#a) and Q~ = P°QP° 1 Q++ = (1 - P°)Q{1 - P°). In 
©, Q is the usual charge operator acting on the Fock space T [551 Eq. (10.52)], 

2 = E a 5(^M/+) - E *o(/r)fc>(/r ), 

i>l i>l 

(fi)i>i and {f~)i>i being respectively orthonormal basis of and . 

Due to (J3J, we have introduced in [25] the notion of P°-trace class operators. We 
say A e & 2 ($)a) is P°-trace class if the operators A++ := (1 - P°)A(1 - P°) and 
A := P°AP° are trace-class (€ @i(^a))j and we define the P°-trace of A by 

(4) trpo A = tr A ++ + tr A . 

Notice, if A is even trace-class then trpo A — tr A. In the following, we denote 
by &i ($)a) the set of all P°-trace class operators. Remark that by definition 

©f(«A) c e 2 (sj A )- 

We have shown in |25l Lemma 2] that any difference of two projectors satisfying 
the Shale-Stinespring criterion, Q = P— P° e &2{$)a), is automatically in 6f ° (5}a)- 
The charge (f2 P |Q|f2p) = trpo(Q) is therefore a well-defined number which indeed 
is always an integer, as proved in |25l Lemma 2]. The P°-trace is an adequate tool 
for describing charge sectors, without using the explicit expression of f2p which can 
be found in the literature. 

In this paper, we study the case of the vacuum: namely we want to show 
the existence of a BDF state fip € B with lowest energy, which we call a BDF- 
stable vacuum. For a small external field, this vacuum will not be charged but 
if the external field is strong enough, we could end up with a charged vacuum, 
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(0p|<2|0p) — tr P o(P — P°) ^ 0. In order to study atoms or molecules, one has to 
minimize the energy in different charge sectors 

B N := {0 P e B | = N} c B. 

In this explained in section 4.2] (see also (23 Remark 6]), the electronic 

orbitals will solve the unprojected Dirac-Fock equations, perturbed by the vacuum 
polarization potentials. It is our goal to study this constrained minimization prob- 
lem in the near future. 

According to Chaix and Iracane Formula (4.1)], the energy of a state Op is 
defined using the renormalized Hamiltonian, acting on the Fock space T , 

(5) M = Jdx: **(x)D«>*(x) : P o +f J ' dx J dy' ^i')^^(v) *■ 

where ^(x) = J2i>i ^{fi)fi{ x )i (fi)i>i being an orthonormal basis of Sj\. The 
choice of the normal ordering with respect to P° corresponds to subtracting the 
energy of the free Dirac sea P° and the interaction potentials involving P°. We 
emphasize that by this choice we make the assumption that the free vacuum is 
unobservable, as done by Dirac UUEj! Heisenberg j23 and Weisskopf [HZj (see 
also (2EJ)- In principle, other choices could be made 

Evaluating the expectation value of Op, we obtain [251 Appendix] 

(6) (0 P |H|Op> = £{Q) 

where Q = P — P° S 6f°(f)A) and £ is the Bogoliubov-Dirac-Fock energy 

(7) £(Q) = tr P a(D°Q) - aD(p Q ,n) + ^D(p Q ,p Q ) - | J J ]Q ^ V ^ dxdy. 
Here Pq{x) = tr^Qix, x) and 

du,,)-*[ 

Jr3 \k\ 2 

Notice that the density pq is well defined due to the ultraviolet cut-off |251 Eq. (9)], 
and that D(f,g) — JJ R6 ^&^H dx dy when / and g are smooth enough. 

As this is seen from (JBJ, the energy of Op only depends on Q = P — P°, which 
is interpreted as the renormalized one-body density matrix of Op. 

3. Existence of a BDF-stable polarized vacuum 

Following a usual method for Hartree-Fock type theories [SSI El El) we may define 
and study the functional £ on the extended convex set 

(8) S A = {Q | < Q + P° < 1, Q e 6f°CS A ), P Q € C}, 

where C is the so-called Coulomb space consisting of functions with finite Coulomb 
norm 

\\p\\l:=D(p,p)=i* ( 

JRZ \ K \ 

More precisely, C is the Fourier inverse of the L 2 space with weight 1/|A:| 2 . 

As our main result we obtain that, for any A, £ is bounded-below and has a 
minimizcr on 5a, therefore there exists a BDF-stable vacuum. 

Theorem 1. Let < a < 4/-7T, n G C. Then £ satisfies, for any Q £ 5a, 

(9) £(Q) + ^D(n,n)>0 
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and it is therefore bounded from below on S\. Moreover, there exists a minimizer 
Q of £ on S\ such that P = Q + P° is a projector satisfying the self-consistent 
equation 

(10) P = X(-oc,o) (D° - atp + ap Q * ± - «^|) ■ 

Additionally, if a and n satisfy 

then this global minimizer Q is unique and the associated polarized vacuum is neu- 
tral: 

(Slp\Q\Slp) =tr P o(Q) =0. 
The proof of this result is given in Section [5j 

Equation (|10|l corresponds to Dirac's picture that the "correct" vacuum P should 
be the projector on the negative spectrum of an effective one-body Hamiltonian. In 
the case without external potential, n = 0, the free projector P° solves l|10f) and is 
the unique BDF-stable vacuum [SJE]- 

Numerically the self-consistent solution of (|10f> could be evaluated by a fixed 
point algorithm, starting with P . In |25j we proved the convergence of this al- 
gorithm to a BDF-stable vacuum solving (|10fl . under reasonable restrictions of the 

form cty/ 1|«||^2 + \\ n \\c — C\ and aVlog A < C2, using the Banach fixed point the- 
orem. This proof is much more constructive than the direct variational proof which 
is given in Section |SJ However, the result of [25] is local in the sense that it is valid 
for weak external potentials ip = n* 1/| • | only. 

The condition Ijlll) means that if the overall charge of the nuclei is not too big 
and a is small enough, the BDF-stable vacuum is unique and stays neutral, cf. 
[221 l^lj . In general, the solution found in Theorem ^ can correspond to a charged 
vacuum. 

There is an interesting symmetry property of the solutions of (|10|) when n is 
replaced by — n. Namely, if P is a solution of i|10|) with external density n, then 
P' = Q' + P° is a solution of l|10|) with external density —n, where Q 1 = —CQC~ X , 
C being the charge conjugation operator |551 page 14]. The two dressed vacua 
P and P' have the same BDF energies and satisfy pq> = —Pq, as suggested by 
the intuition. For this symmetry between matter and antimatter to be true, it is 
essential to have the Fermi level at and not at —1 (see, e.g., the comments of 
page 197] about this fact). 

In Theorem ^ the cut-off A can be chosen arbitrarily large and it is therefore 
natural to describe the behaviour of our solution as A — ^ 00. 

Theorem 2. Let be n e C R L 2 (R 3 ) and < a < 4/n. Then the solution Q\ = 
P\ — P° obtained in Theorem^ satisfies 

|||£T /2 Qa|| 62 ->0, a\\p QA -n\\ c ^0 
as A — > 00, and therefore 

01 

(12) lim minf = D(n,n). 

A^oo Sa 2 

In words, when A — > 00, the vacuum polarization density totally cancels the 
external density n, for pq a — * n in C. But since Qa = — P° — > 0, this means that 
in the limit A — > 00, Qa and pq a become independent. Therefore, the minimization 
without cut-off makes no sense both from a mathematical and physical point of 
view. Indeed Ijl2(l easily implies that when no cut-off is imposed and when ip 7^ 0, 
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the infimum of the functional £ is not attained. In physics, this "nullification" 
of the theory as the cut-off A diverges has been first suggested by Landau et al. 
H3E2HUE1I and later studied by Pomeranchuk et al. (17] . 

In the next section, we propose a renormalization procedure in which we show 
an inequality of the form ^a r logA < 1 where a r is the physical (renormalized) 
coupling constant, different from a. With the usual value a r ~ y^, this leads to 
the physical bound A < 10 280 (in units of mc 2 ). 

The proof of Theorem [21 is given in Section [5] 

Remark. If n is smooth enough, it can be shown that 

III^T /2 0A||e 2 < ^(logA)- 1 , a\\ PQA - n\\ c < ^(logA)- 1 
for some constants C\ and C%. 

4. Reduced energy functional and charge renormalization 

Recall up to now the charge was kept to be the bare one. Next we want to 
derive a renormalization scheme consistent to any order of a for the solution of 
our minimization problem. Note that this procedure is well known in perturbation 
theory, see e.g. ^2 page 194] and 

We first simplify our BDF energy by neglecting the exchange term, 

(13) £ Ted {Q)=tr p0 (D°Q)-a J PqV + ^D(p Q , Pq ). 

From a physical point of view this is quite natural, since the exchange term is usually 
treated together with a term describing the interaction with the photon field to form 
the standard electron self-energy that is a subject of the mass renormalization. 

Notice that since £ re d > £■, the energy functional £ Te d is obviously bounded from 
below on 5a, by Theorem ^ We now state our 

Theorem 3. Let < a < A /it, n € C. Then £ Te d possesses a minimizer Q on Sa, 
which satisfies 

(14) Q = X(-oc.O) (p° ~a(p + ap Q * -L) - P° + 7o , 
where 70 is a finite rank operator of the form 

K 

7o =/ j n i \ip i )((pi\, < rii < 1, 

i=l 

{<Pi)f = i being an orthonormal basis o/ker(_D° — cup + apQ * 1/| • |). 
Additionally, if a and n satisfy 

(15) a^'^^WnWc < 1, 
then this global minimizer Q is unique and 

ker(£>°-a^ + apg*l/|-|)={0} 

which implies 

(16) Q = X(-oo,o)(D -atp + ap Q *y-^ -P°. 

The proof is much simpler than the one of Theorem^ £ rcc i is now a coercive and 
convex continuous functional which is therefore weakly lower semi-continuous on 
the closed convex set 5a, and possesses a minimizer. The proof that it satisfies the 
self-consistent equation l|14H is the same as the one of Theorem except that due 
to the absence of the exchange term, one is not always able to prove that Q + P° 
is a projector, as usual in reduced Hartree-Fock type theories 
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In order to perform our renormalization scheme we expand (|16fl in powers of 
a. Assuming that l|15fl holds, is not in the spectrum of the mean-field operator 
D v + apQ * 1/| • | and we can use the resolvent representation |291 Section VI, 
Lemma 5.6] to derive from 11611 the self-consistent equation for the VP-density 
Pq( x ) = tr C 4Q(x,x) 



(17) PQ(x) = ~[ C 



drj tr C 4 



1 



1 



D° - aip + apQ * pj 



IT) 



IT) 



(x,x). 



1 „ 1 „ 1 



1 „ 1 „ 1 



a'—B—B— + a'—B—B-B- 



1 



Applying the resolvent equation 

1 111 
= a — B — 

A-aB A A A ' ~ A A A A A A A-aB 
and using Furry 's Theorem |17| . telling us that the corresponding a 2 -term with two 
potentials vanish, we obtain 

(18) p Q = aF x [p Q - n] + F 3 [ap Q - an] 

with 



F3[p](x) = 
drj tr C 4 



1 



1 



1 



1 



1 



1 



1 



D° ' | • | D° + " I • I D° + i7i r ' | • | D° - ap + ir]_ 

As realized first by Dirac ^3^2 an( l Heisenberg cf. also JH], the term F\ [p] 
plays a particular role since it is logarithmically ultraviolet divergent. Following, 
e.g., Pauli-Rose |39| . one evaluates in Fourier representation 

F x |>](fc) = -ftk)B A (k), 

with 



-P- 



(x,x). 



(19) 



B A (k) 



z 4 /3 



dz 



1-z 2 l + |/s| 2 (l-z 2 )/4' 

which can be decomposed into ,39. Equ. (5)-(9)] B A (k) — B A — C A (k), with 
(20) 



B A = Sa(0) 
and 



/I+A2 Z 



z 4 /3 



dz 



3tt 



log(A) 



5 

9^ 



3tt 



(21) 



lim C A (k) = C(k) 

A— too 



i r 1 

— / dx(l-x 2 )log[l + fc 2 (l 
2tt Jo 



log2 + 0(l/A 2 ). 



^ 2 )/4], 



which was first calculated by Serber and Uehling • 

Denote p = pq — n the total density, then 118fl reads in terms of p 

(22) p + n = -aB A p + aC A (k)p + F 3 [ap], 
or equivalently 

(23) ap = —an — a 2 B A p + a 2 C A {k)p + aF 3 [ap] 
and 



-F 3 [ap]. 



(24) ap= -— — — n+— — C A (k)ap + — 

1 + aB A 1 + olB a 1 + aB A 

To perform our renormalization scheme we fix as physical (renormalized) objects 
a r p r = ap, with (cf. [13 Equ. (7-18)]) 

(25) = — ^— . 

1 + aB A 

Therefore we can rewrite the self-consistent equation (I22f) as 

(26) a T p x = -a r h + a 2 C A (k)p r + a z F 3 [a T p x ], 
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independently of the bare a. This equation uniquely defines the VP density only 
depending on the physical observable a T , which is what we understand under con- 
sistent to any order. The a r represents the dressed coupling constant, which is 
observed in experiment and whose value is approximately 1/137. Notice that from 
formula Q25JI. it follows that necessarily a r B\ < 1 and a r B A — > 1 as A — > oo. 
We emphasize that although in the literature the expression of a r is sometimes 
expanded to get a r ~ a(l — aB\) leading to the condition aB\ < 1, the real 
constraint indeed applies to the physically observed a r and not the bare one. 

Notice that equation i|26[l satisfied by a r p r is exactly the same as equation (|23|) 
satisfied by ap, but with the logarithmically divergent term a 2 B\p dropped. There- 
fore, as usual in QED |12) . the charge renormalization allows to simply justify the 
dropping of the divergent terms in the self-consistent equation. In practice |37) . one 
would solve (|26[1 with a r ~ 1/137 and with C\(k) replaced by its limit C(k). 



Returning to the effective Hamiltonian D — cap + apQ * 1/ 1 ■ 
i.e. expressing in terms of the physical objects, we obtain 



and inserting l|26|l . 



(27) 



D° + a r p r * — = D° — a T n * r—r + V e g, 



with 



a%CA{k)pr{k) + a r F 3 (a r j9 r 
k 2 



(,•) 



the effective self-consistent potential, where T~ x denotes the inverse Fourier trans- 
form. Notice, this equation is valid for any strength of the external potential. 
However, expanding p r in a T , we obtain to lowest order in a T 



eff 



7T 



a 2 f°° 

H dt(t 



C A (k)h(k) 

k 2 

- _ 1)1/2 



(X) 



2 



i 

t 1 



dx'e' 2 ^- 



t'\t n ( x ') 



the Uehling potential Concerning a point like particle this potential was first 
written down in a closed form by Schwinger |45j . The next term in V e g is of order 
a T (a T Z) 3 . In principle all higher order corrections can be evaluated explicitly, which 
is not the task of our paper. 

Finally we note that the convergence of the term in the right hand side of 1221) , in 
the case of the VP-density in the Furry picture, i.e. aF^an), was shown in various 
papers. The most clarifying proof with respect to spurious third order contributions 
can probably be found in |52| (for earlier proofs, in particular corresponding to 
muonic atoms, we refer to the references in jSJ). However the fact that this term, 
aF^an)^ additionally gives rise to a well defined self-adjoint operator was recently 
proved in |26) . 



5. Proof of Theorems [I] and [21 
In this last Section, we give the proof of our main Theorems. 

5.1. Proof of Theorem ^ The proof that £ is well-defined on Sa can be found 
in details in [23 Theorem 1] . For simplicity, we extend £ to the closed convex set 

S' A = {Q e 6 2 (£a), < Q + P° < 1, p Q e C} 
of the Hilbert space H := {Q & &2(Sja), PQ € C}, by simply letting 
£(Q) - F(Q) - aD(p Q ,n) + ^D(p Q ,p Q ), 
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^^-dxdy ifQe6f°(£ A ) 
6f° 



(28) 

{ +00 ifQ^6f(«A) 

Let us recall the inequality established in 0] 

F(Q) > (l-OTr/4)tr P o(£>°Q) = (1 - a7r/4)(tr(|£»°|Q++) - tv(\D°\Q— )) 

(notice that Q ++ > and Q < when Q £ S' A ), which easily implies the bound 
since 



(29) £(Q) > (l-aTr/4)tr P o(D Q Q) + -\\p Q - 



n \\c ~ ^ll 71 ! 



> ~ n 



This also easily shows that both F and £ are strongly lower semi-continuous and 
coercive on S' A . We now prove that £ is indeed weakly lower semi-continuous (wise) 
on S' A in 7i, which will show the existence of a minimizer since S' A is closed and 
convex, and therefore weakly closed. 

Step 1: £ is wise on S' A . Since the functional 

Q ^ -aD(p Q , n) + -D(p Q , p Q ) = -\\p Q - n\\ 2 c - -\\n\\l 

is easily seen to be wise on the convex set S' A , it only remains to prove that F (defined 
in (|28|l ) is wise on S' A . To this end, we consider a weakly converging sequence Q n — 
Q in 7i, such that Q n € S' A for each n. If liminf„ F(Q n ) = oo, there is nothing 
to show and we can therefore assume that {Qn + )n>i and {Qn~)n>i & re bounded 
in &i($)a)- Due to the cut-off A in Fourier space, (Q n (x, y)) n >i is bounded in the 
Sobolev space 7? 1 (R 6 , C 4 (g) C 4 ), (P|£)0|Q++)n>i and (P|do|q-- )«>i are bounded for 
instance in i? 1 (R 3 , E). We may thus assume, up to a subsequence, that Q n (x, y) — ► 
Q(x,y) in L^ oc (M 6 , C 4 <E> C 4 ), that p| o| Q ++ -> P\d°\q++ and P| D o| Q -- -> Pidg- 
in L 4 oc (K 3 , M) . 

Let us now consider two real functions ry,£ S C°°([0; oo); [0; 1]) such that r/(t) = 1 
if t e [0; 1], ??(*) = if t > 2, < 7y(i) < 1 if t € [1; 2], and + £ 2 = 1. We now 
define tir(x) :— ?7(|a;|/i?) and £r(x) := ^{\x\/R) for a; S M. 3 . In the following, we 
also denote by rjn and £r the multiplication operators by the functions 77^ and £r, 
acting on 5}a- 



Lemma 1. We have II [£r, |D°| 



0(l/R). 



Proof. We compute (?/>| [£,r, |-D°|] |x) in Fourier space, for some tp,x G -Qa (we use 



the notation -E(p) 

(Vlfc R ,|J3 |]|x) = 



U(P ~ q)^P)x(q) (E(q) - E(p)) dpdq 



Els-- 



Els + -))dsdr 



and therefore, using the inequality \E(x) — E{x — y)\ < \y\, we obtain 



\(Mr,\d \}\x)\< 



r£,R(r) dr \\^\\lA\x\\l^ 



and 



< 



»"6?(r) 



dr= R 



r|i(r) 



rZr. 



□ 
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Using this Lemma, we may now write 
tr(|D°|Q++) = tr( V %\D°\Q++)+tr(e R \D°\Q++) 

= tr( m \D \Q+ + VR) + tr(|D |&Q++fr) + tr([fe, \D°\}Q+ + ^ R ) 
= tv{ m \D°\Q++ m ) + tr(|D°|^Q++^) + 0(1/R) 



since 



\tr(^ R ,\D°\}Q+ + U)\ < ||[^,|-D |]|| Soo(i5A) ||Q+ + || ei(AA) =0(l/ii:) 

by Lemma ^ and since by assumption (Q^ + ) n >i is bounded in &i(Sj\). With the 
same argument for Q~~ , we obtain 

tr P o(£>°Q„) = tr( m \D°\Q++ VR ) - tr( m \D°\Q-~ m ) 

+ tr(\D°\Z R Q++Z R ) - tr(\D°\Z R Q--Z R ) + 0(1/ R). 

On the other hand, we have 



\x-y\ J J k — 2/| 

2ln t„ „,M2 



JJ \x-y\ 



\x - y\ 

since 

m(x) 2 ^ R (y) 2 \Q n (x,y)\ 2 H^III^a) 



\x-y\ dxdy ~ R ' 

We therefore obtain 

F(Q n ) = ti(r, R \D°\Q++r, R ) - ix(ti R \lft\Q~ri R ) 

a f f VR(x) 2 V3R(y) 2 \Qn(x 1 y)\ 2 , , . , nn o, t s ■ /m 0| t r, — t n 

>' : dxdy + tr(\D v \£ R Q^£ R ) - tr(|Z7 \£ R Q n £ R ) 

f - y\ 

a f f U(x) 2 \Q n (x,y)\ 2 



J J \x-y\ 



Notice now that < Q n + P° < 1 implies |Q„| 2 < Q++ - Q~~ (see 0J). We now 
localize this inequality to obtain £, R \Q n \ 2 £,R < £,RQt + £,R — £,RQn ^R- By Kato's 
inequality [1], we now have 

^R{x) 2 \Q„(x,y)\ 2 7T , ln0|e 1 , 

cfecfy < — tr(|_D |G?.<3„60 



|a: -I/I - 2 

7T 

2 



< J(tr(|i? |^Q+ + 60 " tr(\D°\Z R Q--{[ R )) 



and therefore, since (1 — air /4) > 0, 

F(Q„) > tr(r ?fl |£>°|Q++» ? R) - tr(» J H|D |Q--» JJi ) 

_ a re m(x) 2 mR {y) 2 \Q n { X ,y)\ 2 ^ + 
2JJ \x-y\ 
Passing now to the limit as n — > oo and using the local compactness of Q n (x. y) in 
L 2 oc (R 6 ) and P\ D o lQ ++, P\d"\q-- in L ioc( R3 )> we obtain 

liminf F(Q n ) > tt(r) R \D°\Q ++ VR ) - tr(r, R \D°\Q— r, R ) 

-2JJ \x~^\ dxdy + 0(1/ R). 

If we now let R — > oo, we obtain liminf „^oo F(Q n ) > F(Q) and therefore F is wise 
on S' A . 
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Step 2: at least one of the minimizers satisfies <|1(J[I . In the previous step, we 
have shown the existence of a minimizer. It now remains to show that one of them 
indeed satisfies <|10|) . 

Lemma 2. Let Q be a minimizer of £ in S' A . Then either Q + P° is a projector, 
or 

(30) Q + P° = P + fj,\f)(f\, 

where P is a projector, /i € (0; 1) and f G kei(DQ), with 

Dq:= D -aip + a P Q * — - - a- r. 

v v hi \x-y\ 

Proof. Our proof is inspired by classical arguments already used in the Hartree-Fock 
theory EJ. 

Notice that since Q is compact, Q + P° is a compact perturbation of P° and 
therefore its essential spectrum is cr C ss(0 + P°) = {0, 1}, meaning that a(Q + P°) fl 
(0; 1) only contains eigenvalues of finite multiplicity accumulating at {0, 1}. Let us 
assume that Q + P° possesses two different eigenvectors ipi , ip 2 ■ 

Q + P° = \M(<pi\ + \ 2 \<P2)(<P2\ + G 

where Ai, A2 <G (0; 1) and Gipi = Gif2 = 0. We now introduce Q t := Q + e\(fi) (<pi\ — 
e| V2) (^2 1 which belongs to 5a for e small enough and compute 

£(Q e ) = £(Q) + e((^\D Q \^) - (tp 2 \D Q \<p 2 )) ~ <?~ j j l -^0^dxdy. 

Therefore, using either the first order term in e if it does not vanish, or the second 
order term, we can always decrease the energy. This is a contradiction which implies 
that a(Q + P°) n (0; 1) contains at most one eigenvalue of multiplicity 1 and thus 

Q + P° = P + »\f)(f\ 

where P is a projector and /i e [0; 1). If /i =/= 0, using the same type of variation 
Q e := Q + e|/)(/|, we easily show that indeed / € ker(Z?Q). □ 

If Q is a minimizer of the form (|30[) . we now see that 

£{Q) = £(P-P°)+fi(D P _ P ofJ) 
= £(P-P°) + v(D Q f,f) 
= £(P-P°) 

and therefore P — P° is also a minimizer of £ (i.e. P is BDF-stable vacuum). In 
|25l proof of Theorem 2] , we have already shown that a minimizer of £ on 5a taking 
the form P — P° where P is an orthogonal projector, is indeed a solution of the 
self-consistent equation (|10() . 

Step 3: uniqueness of the global minimizer of £ under the condition . Due to 
[251 Theorem 2] , we know that the global minimizer Q of £ is unique if Dq satisfies 

d\D Q \ > \D°\ 

for some d such that adn/A < 1. 

We know that £ (Q) < £ (0) =0 and therefore, by an argument similar to (|29|l . 

(3D ( 2 --fj II l -^fd X dy + f||p, -n||g<f Ml 

^JJJsfi \x-y\ 2 2 

and thus 

(32) // ^ 2 dxd y< ^ — 



1 - an/4 
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(33) \\pQ-n\\c<\\n\\ c . 

Recall that Dq = D° + atp'^ — aRg where <p'q — (pq — n) * 1/| • | and Rq is the 
operator with kernel Q(x, y)/\x — y\. Now, we have 



< 



<Pq 



\D°\ 



<(2tt) 



-1/2 



6 6 (f> A ) 



\ E (-r l \ 



where we recall that E(p) = yl + p 2 , and by Theorem 4.1]. Therefore 
, 1 



\D°\ 



<S 6 



< SqCq 



(47r)5 6 C 6 ||p -n|| c 



with S 6 = 2- 5 / 6 3 1 / 6 tt- 1 /6 and where Cf . = 3-1/622/3^-2/3 is the Sobolev constant 
for the inequality ||/||l 6 (k 3 ) < Cell V/||l2( K 3). Due to this shows that 

|^|<«IH| C \D°\, 

where k = (47r)S 6 C 6 = 7r 1/6 2 11/6 . On the other hand, we know from [211 Proof of 
Lemma 4] that 



Rn < \ ~ 



7T ff \Q(x,yW 



F - 2/1 



and therefore, using (|52l . 



I^qI < - 



a/2 



2 V 1 - an /4 



Inllc |D°I 



As a conclusion, when 



a/2 



2 V 1 - a7r/4 



k \\n\\ c < 1, 



Dq fulfills d|£»g| > |D°| with 



d= <1- a - 



a/2 



2 y 1 - a?r/4 



+ K ||n||c 



Applying now |25l Theorem 2], we obtain that the minimizer Q is unique when 
adir/4 < 1, i.e. under the condition (HJ. 

Assuming now that holds, let us show that the unique BDF- stable vacuum 
P is not charged. To this end, we define, for t £ [0; 1], 



Q(t) = X(-oo-o) D" + at{p Q - n) * 



at 



Q{x,y) 
\x - y\ 



t t— ► Q(t) is a continuous function for the &2(^)a) topology, since by the previous 
estimates Dq(t) = D° + at(pQ — n) * A — at ®x-y\ possesses a gap around 0, 
uniformly in t € [0; 1]. This implies that 

q :t^tT P o{Q(t)) = tr(Q(t) 3 ) 

is continuous on [0; 1], by |251 Lemma 2]. Since q(Q) — and q(t) is an integer for 
any t £ [0; 1], we therefore deduce that 

q(l) = (Slp\Q\(lp) =trpo(Q) = 0. 

This ends the proof of Theorem ^ D 
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5.2. Proof of Theorem |5J We first prove l|12fl which will easily imply that Q A 
obtained by Theorem ^ behaves at stated as A — > oo, due to l|29|) . To this end, we 
introduce 

D° - an A * -i- ) - P° € 5 A , 



!A ■ = X(-oo;0) 

n(k) 



1 + aB A {k) 

where we recall that B A (k) = B A — C A (k) is defined in l|19|) . We now show that 

a 

lim £ (Q A ) = --D(n,n), 

A— >oo z 

which will imply (fT^j) , by 10 . 

Let us now compute p& '•= Pq a - By Pa satisfies 

^ A ( fc ) = T~ p ^n ^W^aO) +-F3 a«A (fe), 

and therefore 

pA(A;) - = -n\(k) + F 3 [an A ](k). 
When a > 0, since (1 + qBa(/c)) _1 — > a.e., we obtain by Lebesgue's dominated 
convergence Theorem that ||n-A||cnL 2 — > as A — > 00. By the fixed-point estimates 
of in the case of the reduced model (they are then independent on the cut-off 
A as this can be seen from the proof of |251 Theorem 3]), it is known that F3 is 
continuous at for the C C\ L 2 topology. We therefore obtain 

lim ||pa - n\\ CnL 2 = 0. 

A^oo 

On the other hand, we also know from the bounds proved in |25| . that 

trpo^QA) 1 ^ = tv(\D \Ql)^ < Ca\\ PA - n\\ cnL , 
for some constant C independent of A. Therefore 

(X 

lim £(Qa) = -—D(n,n) 

A— too Z 

which ends the proof of Theorem □ 

References 

1. A. A. Abrikosov, I. Halatnikov and L.D. Landau. On the quantum theory of fields. Nuovo Cim. 
Suppl, 3: 80-104, 1956. Reprinted in Collected papers of L.D. Landau, edited by D. Ter Haar, 
Pergamon Press, 1965. 

2. A. Aste, G. Baur, K. Hencken, D. Trautmann and G. Scharf. Electron-positron pair production 
in the external electromagnetic field of colliding rclativistic heavy ions. Eur. Phys. J. C, 23(3): 
545-550, 2002. 

3. V. Bach. Error bound for the Hartree-Fock energy of atoms and molecules. Comm. Math. 
Phys., 147: 527-548, 1992. 

4. V. Bach, J.-M. Barbaroux, B. Helffer and H. Siedentop. On the Stability of the Relativistic 
Electron-Positron Field. Commun. Math. Phys., 201: 445-460, 1999. 

5. V. Bach, E.H. Lieb and J. P. Solovej. Generalized Hartree-Fock theory and the Hubbard model. 
J. Statist. Phys., 76 (1-2): 3-89, 1994. 

6. E. Borie and G.A. Rinker. Rev. Mod. Phys, 54: 67 (1982) 

7. P. Chaix and D. Iracane. From quantum electrodynamics to mean field theory: I. The 
Bogoliubov-Dirac-Fock formalism. J. Phys. B. , 22(23): 3791-3814, 1989. 

8. P. Chaix, D. Iracane and P.L. Lions. From quantum electrodynamics to mean field theory: II. 
Variational stability of the vacuum of quantum electrodynamics in the mean-field approxima- 
tion. J. Phys. B., 22(23): 3815-3828, 1989. 

9. P. Chaix. Une Methode de Champ Moyen Relativiste et Application a I'Etude du Vide de 
I'Electrodynamique Quantique. PhD Thesis, University Paris VI, 1990. 

10. P.-A.-M. Dirac. Theorie du positron. Solvay report, pages 203—212. Paris: Gauthier-Villars. 
XXV, 353 S., 1934 (reprinted in Selected papers on Quantum Electrodynamics, edited by J. 
Schwinger, Dover, 1958). 

11. P.A.M. Dirac. Discussion of the infinite distribution of electrons in the theory of the positron. 
Proc. Camb. Philos. Soc, 30: 150-163, 1934. 



SELF-CONSISTENT POLARIZED VACUUM 



15 



12. F. J. Dyson. The S Matrix in Quantum Electrodynamics. Phys. Rev., 75 (11): 1736-1755, 
1949. 

13. M.J. Esteban, E. Sere. Solutions of the Dirac-Fock Equations for Atoms and Molecules. Comm. 
Math. Phys., 203: 499-530, 1999. 

14. M.J. Esteban, E. Sere. Nonrclativistic limit of the Dirac-Fock equations. Ann. Henri Poincare, 
2(5): 941-961, 2001. 

15. M.J. Esteban, E. Sere. A max-min principle for the ground state of the Dirac-Fock functional. 
Contemp. Math., 307: 135-139, 2002. 

16. L.L. Foldy and E. Eriksen. Some Physical Consequences of Vacuum Polarization. Phys. Rev., 
95(4): 1048-1051, 1954. 

17. W.H. Furry. A symmetry theorem in the positron theory. Phys. Rev., 51: 125, 1937. 

18. W.H. Furry. On Bound States and Scattering in Positron Theory. Phys. Rev., 81(1): 115-124, 
1951. 

19. W.H. Furry and J.R. Oppenheimer. On the Theory of the Electron and Positive. Phys. Rev., 
II. Ser., 45: 245-262, 1934. 

20. H. Ficrz and G. Scharf. Particle interpretation for external field problems in QED. Helv. Phys. 
Acta, 52: 437-453, 1979. 

21. R. Glauber, W. Rarita, and P. Schwed. Vacuum polarization effects on energy levels in fi- 
mesonic atoms. Physical Review, 120(2): 609—613, 1960. 

22. W. Greiner, B. Miiller, and J. Rafclski. Quantum Electrodynamics of Strong Fields. Texts and 
Mongraphs in Physics. Springer- Verlag, lrst edition, 1985. 

23. A. Hamm and D. Schiitte. How to remove divergences from the QED-Hartrce approximation. 
J. Phys. A: Math. Gen., 23: 3969-3982, 1990. 

24. C. Hainzl. On the Vacuum Polarization Density caused by an External Field. Ann. Henri 
Poincar, 5(6): 1137-1157, 2004. 

25. C. Hainzl, M. Lewin and E. Sere. Existence of a stable polarized vacuum in the Bogoliubov- 
Dirac-Fock approximation, to appear in Commun. Math. Phys. 

26. C. Hainzl and H. Siedentop. Non-Pcrturbative Mass and Charge Rcnormalization in Rclativis- 
tic no-photon Quantum Electrodynamics. Comm. Math. Phys., 243: 241-260, 2003. 

27. W. Heisenberg. Bemerkungcn zur Diracschcn Theorie des Positrons. Z. Phys., 90: 209-231, 
1934. 

28. C. Itzykson and J.-B. Zuber. Quantum Field Theory. McGraw-Hill, New York, 1980. 

29. T. Kato. Perturbation Theory for Linear Operators, volume 132 of Grundlehren der mathe- 
matischen Wissenschaften. Springer- Verlag, Berlin, 1 edition, 1966. 

30. M. Klaus and G. Scharf. The regular external field problem in quantum electrodynamics. Hev. 
Phys. Acta, 50: 779-802, 1977. 

31. M. Klaus and G. Scharf. Vacuum Polarization in Fock Space. Helv. Phys. Acta, 50: 803-814, 
1977. 

32. L.D. Landau. On the quantum theory of fields. Bohr Volume, Pergamon Press, Oxford, 1955. 
Reprinted in Collected papers of L.D. Landau, edited by D. Ter Haar, Pergamon Press, 1965. 

33. L.D. Landau and I. Pomcranchuk. On point interactions in quantum electro dynamics. Dokl. 
Akad. Nauk SSSR, 102: 489-492, 1955. Reprinted in Collected papers of L.D. Landau, edited 
by D. Ter Haar, Pergamon Press, 1965. 

34. L.D. Landau. Fundamental problems. Pauli Memorial Volume, Interscience, p 245, 1960. 
Reprinted in Collected papers of L.D. Landau, edited by D. Ter Haar, Pergamon Press, 1965. 

35. E.H. Lieb. Variational Principle for Many-Fermion Systems. Phys. Rev. Lett., 46: 457-459, 
1981. 

36. E.H. Lieb and H. Siedentop. Renormalization of the regularized relativistic electron-positron 
field. Comm. Math. Phys., 213 (3): 673-683, 2000. 

37. P. J. Mohr, G. Plunicn and G. Soff. QED corrections in Heavy Atoms. Phys. Rep, 293: 227- 
369, 1998. 

38. E. Paturcl. Solutions of the Dirac-Fock equations without projector. Ann. Henri Poincare, 
1(6): 1123-1157, 2000. 

39. W. Pauli and M.E. Rose. Remarks on the Polarization Effects in the Positron Theory. Phys. 
Rev II, 49: 462-465, 1936. 

40. I.Ya Pomeranchuk, V.V. Sudakov and K.A. Ter-Martirosyan. Vanishing of rcnormalizcd 
charges in field theories with point interaction. Phys. Rev., 103 (3): 784-802, 1956. 

41. P.-G. Reinhard, W. Greiner and H. Arenhvel. Electrons in Strong External Fields. Nucl. Phys. 
A, 166: 173-197, 1971. 

42. J. Rcinhardt, B. Miiller and W. Greiner. Theory of positron production in heavy-ion collision. 
Phys. Rev. A, 24(1): 103-128, 1981. 

43. S.N.M. Ruijsenaars. On Bogoliubov transformations for systems of relativistic charged parti- 
cles. J. Math. Phys., 18(3): 517-526, 1977. 



16 



C. HAINZL, M. LEWIN, AND E. SERE 



44. G. Scharf and H.P. Scipp. Charged Vacuum, Spontaneous Positron Production and All That. 
Phys. Lett., 108B(3): 196-198, 1982. 

45. J. Schwinger. Quantum Electrodynamics II. Vacuum Polarization and Self-Energy. Phys. Rev., 
II. Ser., 75: 651-679, 1949. 

46. J. Schwinger. On the Green's Function of Quantized Fields. II. Proc. Nat. Acad. Sci., 
37(7):455-459, 1951. 

47. H.P. Seipp. On the S-operator for the external field problem of QED. Helv. Phys. Acta, 55: 
1-28, 1982. 

48. R. Serber. Linear modifications in the Maxwell field equations. Phys. Rev., II. Ser., 48: 49-54, 
1935. 

49. V.M. Shabaev. Two-time Green's function method in quantum electrodynamics of high-Z 
few-electron atoms. Phys. Rep., 356, 119-228 (2002) 

50. D. Shale and W.F. Stincspring. Spinor representations of infinite orthogonal groups. J. Math, 
and Mech., 14: 315-322, 1965. 

51. B. Simon. Trace Ideals and their Applications. Vol 35 of London Mathematical Society Lecture 
Notes Series. Cambridge University Press, 1979. 

52. G. Soff and P. Mohr. Phys. Rev. A, 38, 5066 (1988) 

53. J. P. Solovej. Proof of the ionization conjecture in a reduced Hartrcc-Fock model. Invent. Math., 
104(2): 291-311, 1994. 

54. B. Swirles. The relativistic self-consistent field. Proc. Roy. Soc. A, 152: 625-649, 1935. 

55. B. Thaller. The Dirac Equation. Springer Verlag, 1992. 

56. E.A. Uchling. Polarization effects in the positron theory. Phys. Rev., II. Ser., 48: 55-63, 1935. 

57. V. Weisskopf. Uber die Elcktrodynamik des Vakuums auf Grund der Quantentheorie des Elek- 
trons. Math.-Fys. Medd., Danske Vid. Selsk., 16(6):l-39, 1936 (reprinted in Selected papers 
on Quantum Electrodynamics, edited by J. Schwinger, Dover, 1958). 

CEREMADE, Universite Paris-Dauphine, Place du Marechal de Lattre de Tassigny, 
F-75775 Paris Cedex 16, France & Laboratoire de Mathematiques Paris-Sud-Bat 425, 
F-91405 Orsay Cedex 

E-mail address: hainziaceremade.dauphine.fr 

CEREMADE, Universite Paris-Dauphine, Place du Marechal de Lattre de Tassigny, 
F-75775 Paris Cedex 16, France. 

E-mail address: lewinSceremade.dauphine.fr 

CEREMADE, Universite Paris-Dauphine, Place du Marechal de Lattre de Tassigny, 
F-75775 Paris Cedex 16, France. 

E-mail address: sereaceremade.dauphine.fr 



